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ABSTRACT 

Pair-separation statistics of in-situ and synthetic surface drifters deployed near the Deepwater 
Horizon site in the Gulf of Mexico are investigated. The synthetic trajectories derive from a 1-km- 
resolution data-assimilative Navy Coastal Ocean Model (NCOM) simulation. The in-situ drifters 
were launched in the Grand LAgrangian Deployment (GLAD). Diverse measures of the dispersion 
are calculated and compared to theoretical predictions. For the NCOM pairs, the measures indi¬ 
cate nonlocal pair dispersion (in which pair separations grow exponentially in time) at the smallest 
sampled scales. At separations exceeding 100 km, pair motion is uncorrelated, indicating absolute 
rather than relative dispersion. With the GLAD drifters however the statistics are ambiguous, with 
some indicating local dispersion (in which pair separations exhibit power law growth) and others 
suggesting nonlocal dispersion. The difference between the two data sets stems in part from inertial 
oscillations, which affect the energy levels at small scales without greatly altering pair dispersion. 

These were significant in GLAD but much weaker in the NCOM simulation. In addition the GLAD 
drifters were launched over a limited geographical area, producing few independent realizations and 
hence lower statistical significance. Restricting the NCOM set to pairs launched at the same loca¬ 
tions yields very similar results, suggesting the model is for the most part capturing the observed 
dispersion. 


1. Introduction 


Submesoscale processes, i.e., with length scales of 


0.1-10 km (Thomas et al. 2008), are believed to be 


important in the upper ocean (McWilliams 2008 
Klein 2009). These are the transition scales be¬ 


tween the largely balanced quasi-2D flows at the 
mesoscales and 3D (unbalanced) flows at smaller 
scales. While observational evidence of submesoscale 
activity in the ocean is accumulating, important 
questions about their dynamics and the consequences 
for transport remain. 

In the atmosphere, the balanced scales are char¬ 
acterized by a kinetic energy spectrum proportio nal 
to k~ 3 , where k is the horizontal wavenumber (Nas- 
trom and Gage |T985 ). The consensus is that this 
reflects a quasi-2D enstrophy cascade toward smaller 
scales (Kraichnan |1967 Charney1|1971 ). At scales of 
several hundred kilometers, the spectrum transitions 
to a A; -5 / 3 dependence. This primarily reflects diver- 
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gent motions (inertia-gravity waves), at scales where 
the Rossby number exceeds one ( Callies et al.|2014 ). 
Callies and Ferrari (|2013 ) suggest a similar situation 


exists in the ocean. 

The slope of the energy spectrum is important for 
Lagrangian transport (e.g., Bennett| 2006 LaCasce 


20081. With a —5/3 slope the dispersion of pairs of 
particles (or “relative dispersion”) is local, meaning 
separations between pairs of particles are dominated 
by eddies of comparable scales. With a —3 or steeper 
slope the dispersion is nonlocal and governed by the 
largest eddies in the k~ 3 range. Local dispersion re¬ 
sults in small scale “billowing,” as with smoke from a 
stack, while nonlocal dispersion produces filaments. 
Particle dispersion can thus be used to infer aspects 
of the energy spectrum, which can be useful in the 
ocean at scales below those resolved by satellite al¬ 
timetry. 

Hereafter we examine relative dispersion at the 
surface of the Gulf of Mexico (GoM). The study 
was motivated by the Grand LAgrangian Deploy¬ 
ment (GLAD), which was conducted in the vicinity 
of the Deepwater Horizon (DwH) site in July 2012 
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and in which a large number of surface drifters were 


deployed (Olascoaga et al. 

2013 

Poje et al. 

2014 

Jacobs et al. 

2014 

Coelho et al. 

2015 

). A primary 


goal of GLAD was to study dispersion at the subme- 
soscales in the GoM. 

Relative dispersion in the GoM has been stud¬ 


ied previously. LaCasce and Ohlmann (2003) exam¬ 


ined “chance pairs” of drifters (i.e., drifters not de¬ 
ployed together) from the Surface-CUrrent and La- 


grangian drifter Program (SCULP) (Ohlmann and 
Niiler|2005 1 and found nonlocal dispersion below the 
deformation radius, Ld, which is approximately 45 
km in the GoM (Chelton et al. 19981. Support¬ 


ing evidence, using pair separation probability distri¬ 


bution functions (PDFs), was obtained by LaCasce 


(2010). However, using different measures (the sec¬ 


ond order longitudinal velocity structure function 
and the separation-averaged relative diffusivity) with 


the GLAD drifters, Poje et al. (2014) concluded the 


dispersion was local, from few hundred meters to sev¬ 
eral hundred kilometers, implying a shallower kinetic 
energy spectrum. 

Relative dispersion is often studied using two types 
of measures (LaCasce 20081. The first treats time 


as the independent variable. This includes the rel¬ 
ative dispersion (the mean square pair separation), 
the kurtosis (the normalized fourth moment) and the 
separation PDF. The second type uses distance as 
the independent variable. This includes the structure 
functions, the separation-averaged relative diffusivity 


and the finite-scale Lyapunov exponent (FSLE) (Ar- 


tale et al. 1997 Aurcll et al. 1997). The two types 


often produce different results, but the reasons for 
this are rarely examined. 

In the present paper, we examine both types of 
statistics, using synthetic particles and drifters in the 
GoM. The former were obtained by integrating sur¬ 
face velocities produced by a data-assimilative simu¬ 
lation with the 1-km-resolution Navy Coastal Ocean 


Model (NCOM) (Jacobs et al. 20141, and the lat¬ 
ter are the drifters in the GLAD experiment. The 
model allows for large numbers of particles, increas¬ 
ing statistical reliability, whereas the drifters more 
accurately reflect the actual situation in the GoM. 

The paper is organized as follows. In Section 2 
and Appendix A we present relevant theory for pair- 
separation statistics. In Section 3 we examine the 
NCOM pair separations, and the GLAD pairs in Sec¬ 
tion 4. A summary and concluding remarks are of¬ 
fered in Section 5. Details of the numerical simula¬ 
tion and the GLAD experiment are given in Appen¬ 
dices B and C, respectively. 


2. Theory 

Let ro be the distance between two fluid parcels 
at time t = 0 and r be the separation at time t. 
In homogeneous, stationary, and isotropic 2D turbu¬ 
lence, the PDF of pair separations, p(r,t ), obeys a 
Fokker-Planck equation: 

d t p = r~ 1 dr{rn 2 drp ), (1) 


where k 2 (r) is the scale (r) dependent relative dif¬ 
fusivity. The 3D version of 0 was proposed by 


Richardson (1926) to describe smoke dispersion in 
the atmospheric boundary layer. Richardson ob¬ 
tained a self-similar solution, based on an empiri¬ 
cal diffusivity derived from observations. The equa¬ 
tion was later derived by Kraichnan (1966) using his 


“abridged Lagrangian history direct interaction ap¬ 
proximation” and by Lundgren (1981), assuming an 


advecting velocity with a short correlation time. For 

Chapter 11). 


an overview, see 


Bennett (2006 


Pair dispersion depends on whether the pair veloc¬ 
ities are correlated or not, specifically whether the 
normalized Lagrangian velocity correlation 


2 (vi ■ Vj ) 
<«?> + <«?> 


= 1 - 


({Vj - Vj ) 2 ) 

<«?> + <«?> 


( 2 ) 


(where the angle bracket indicates statistical aver¬ 
age) equals 1 or 0, respectively (e.g., 

2009). The second term on the right 
tional to the second-order velocity structure func¬ 
tion, 

5 2 (r):=(u 2 )^((u l -u J ) 2 ), (3) 

where v is the difference in the Eulerian velocity 
between points separated by a distance r. The 
Eulerian-Lagrangian equivalence in (|3]) is a distin¬ 
guishing aspect of homogeneous, isotropic turbulence 
(B ennett||1984 ). 

At large separations, when the pair velocities are 
uncorrelated, the relative diffusivity k 2 is constant 
and equal to twice the single particle diffusivity. At 
smaller scales, n 2 can be inferred if the energy spec¬ 
trum has a power law dependence, i.e., E(k) oc k~ a 


(|Bennett| 

1984 

). In either case equation (fl|) can then 

be solved 

(Bennett 2006 

). Relevant 2D solutions are 

given in I 

^aCasce 

(2010 

and 

Graff et al. (2015), and 


are reproduced in Appendix A. From these, the (raw) 
statistical moments, given by 


rOO 

(r n ) := 2 tt / r n+1 p(r, t) dr, (4) 

Jo 

can be calculated. Table [T] shows the time de¬ 
pendences for (r 2 ) (the relative dispersion) and the 
(r 4 )/(r 2 ) 2 (kurtosis) in the three specific 2D disper¬ 
sion regimes considered here. 


Koszalka et al. 
side is propor- 
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Table 1. The pair-separation measures in the three regimes 
considered here. The measures are the kinetic energy spec¬ 
trum, E(k ), where k is wavenumber; the two particle diffu- 
sivity, K 2 (r), where r denotes separation; the second-order 
longitudinal structure function, 52 M; the relative disper¬ 
sion, ( r 2 ); and the separation kurtosis, (r 4 )/(r 2 ) 2 . The latter 
two are based on the (raw) statistical moments of the time- 
dependent probability distribution function (PDF) of pair sep¬ 
arations, obeying 0 . The results shown for the Richardson 
and Rayleigh cases are the asymptotic (long time) limits. 



Lundgren 

Richardson 

Rayleigh 

E{k) 

oc k~ 3 

OC fc -5 / 3 

- 

K2 ( r) 

f~y 

<N 

s- 

II 

= dr 4 / 3 

= const 

S 2 (r) 

oc r 2 

oc r 2 / 3 

= const 

<r 2 > 

— r 2 e 8 t/T 

~ 5.2675 /3 3 t 3 

~ 4/^2^ 

<r 4 )/(r 2 ) 2 

= e 8 ‘/ T 

~ 5.6 

~ 2 


The nonlocal regime [which we refer to as the 
“Lundgren regime,” after Lundgren (1981)] corre¬ 
sponds to an energy spectrum at least as steep as 
k~ 3 , with a structure function S 2 oc r 2 . The PDF is 
not self-similar, but rather becomes more and more 
peaked. Both the dispersion and kurtosis increase ex¬ 
ponentially in time. The local regime we will consider 
has a Kolmorogorov energy spectrum, E oc fc -5 / 3 , or 
equivalently S 2 oc r 2 / 3 , and referred as the “Richard¬ 
son regime” after Richardson (1926). In this case the 
separation PDF asymptotes to a self-similar form, 
with a kurtosis of 5.6, and the dispersion increases 
as time cubed. With uncorrelated pair velocities, 
the second-order structure function is constant with 
separation (and equal to twice the mean square sin¬ 
gle particle velocity). The PDF also asymptotes to 
a self-similar form, with kurtosis of 2 and the dis¬ 
persion increasing linearly in time. This self-similar 
PDF is a Rayleigh distribution, so we refer to this as 
the “Rayleigh regime.” 


3. Simulated pair-separation statistics 

The simulated trajectories were constructed by in¬ 
tegrating surface velocities produced by an NCOM 
simulation (cf. Appendix B). The integrations were 
carried out using a stepsize-adapting fourth/fifth- 
order Runge-Kutta method with interpolations ob¬ 
tained using a cubic scheme. One-month-long 
records, with 10 positions per day, were produced 
with a range of initial separations, from the smallest 
scale resolved by the model simulation up to 30 km 
(ro = 1, 5, 10 and 30 km). The trajectories were ini¬ 
tiated every other day in the northern GoM near the 
DwH site, in two 5x5 100-km-width grids displaced 


30.0°N 


28.8°N 


27.5°N 


26.3°N 


25.0°N 

90.0°W 88.8°W 87.5°W 86.2°W 85.0°W 

Fig. 1. The rectangles indicate the boundaries of the grids 
of initial positions for the integration of trajectory pairs (with 
ro = 1 km initial separation) using surface velocities produced 
by the 1-km-resolution Navy Coastal Ocean Model (NCOM) 
simulation of the Gulf of Mexico (GoM). Dots indicate deploy¬ 
ment locations of drifters from the Grand LAgrangian Deploy¬ 
ment (GLAD). 



by the chosen separation. The reference and auxil¬ 
iary grids with ro = 1 km are shown in Fig. [lj The 
trajectories were started 1 July 2013 and 1 February 
2014 to survey summer and winter conditions. The 
two were expected to exhibit different pair-separation 
statistics, as the mixed layer is deeper in winter. 

However, snapshots of the instantaneous surface 
vorticity (Fig. [2]) reveal roughly the same range of 
eddy scales in the two seasons. The results seen 
hereafter similarly show only small changes with sea¬ 
son. These figures also reveal that the West Florida 
Shelf and the Bay of Campeche are relatively eddy 
inactive regions, and that an anticyclonic ring has 
pinched off from the Loop Current. Apart from these 
regions however the eddy field does not exhibit sig¬ 
nificant spatial variability. Thus we assume homo¬ 
geneity holds fairly well, as assumed in Section [2j 

An additional theoretical prerequisite is station- 
arity. While temporal variability of the simulated 
background eddy field is evident, this mainly mani¬ 
fests on seasonal timescales. Therefore, considering 
motion over a period of 1 month, as we do here, is 
not restrictive but rather necessary for stationarity 
to be fairly well guaranteed. 

Isotropy, the remaining prerequisite, is also real¬ 
ized. This can be seen by plotting the ratio of the 
zonal to the meridional relative dispersion ratio as 
a function of scale ( |Morel an d Larcheveque 1974), 
as shown in Fig. [3] Isotropy holds in both summer 
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Fig. 2. Snapshots of surface vorticity (normalized by the 
mean Coriolis parameter in the GoM) from the NCOM simu¬ 
lation in summer (top) and winter (bottom). 



r [km] 


Fig. 3. Zonal-to-meridional relative dispersion ratio as 
a function of pair separation in summer (solid) and winter 
(dashed) based on simulated trajectory pairs with various ini¬ 
tial separations. The vertical line indicates the gravest baro- 
clinic Rossby deformation radius. 


(solid) and winter (dashed), irrespective of the ini¬ 
tial pair separation. At separations exceeding 100 
km, the dispersion becomes zonally anisotropic, as 
it does in the atmosphere (Graff et al. 20151. But 
below that, the dispersion is isotropic. 

We then determine over which scales the pair mo¬ 
tion is correlated. Inspection of Fig. [4] reveals that, 
fairly insensitive to the initial pair separation, the 
motion is correlated below the deformation radius, 



Fig. 4. The normalized Lagrangian velocity correlation as 
function of scale in summer (solid) and winter (dashed) for 
various initial separations in the NCOM simulation. The ver¬ 
tical line indicates the Rossby deformation radius. 


in both seasons. At Ld, the correlation is roughly 
0.5. The motion is strongly correlated below 20 km 
and decorrelated above 100 km. Thus the proper 
framework for interpretation above 100 km is abso¬ 
lute dispersion, possibly with the inclusion of a mean 
zonal shear. 

We now inspect the pair separation PDFs and 
compare them to the theoretical predictions (Fig. [5]). 
We calculated the PDFs using kernel density estima¬ 
tion as implemented in Matlab’s function ksdensity 
with positive support. Each theoretical PDF has two 
parameters: the initial separation, r 0 , and a growth 
parameter. Since all pairs have the same initial sep¬ 
aration, we assume tq is the same as in the sim¬ 
ulation. The growth parameters for the Lundgren 
and Richardson distributions (T and /3, respectively) 
were determined by fitting the dispersion using least 
squares during the initial period, up to the point 
when the root-mean-squared (rms) separation was 
a factor a greater than the initial value. We chose 
a = 5, but the results are fairly insensitive to the 
choice. The parameter for the Rayleigh PDF (k 2 ) 
was calculated from the relative dispersion at late 
times, after the pair motion was decorrelated, as de¬ 
scribed in Appendix A. 

The PDFs are plotted with the theoretical curves 
at t = 3.5 d. Several initial separations are used, 
for particles deployed in summer (in red) and win¬ 
ter (in blue), and the parameters obtained for the 
theoretical curves are indicated in the inserts. The 
observed and theoretical PDFs were compared us¬ 
ing the Kolmogorov-Smirnov (KS) statistic (cf., e.g., 
; since the winter and summer dis¬ 
similar, we used the summer PDF 
for this. The degrees of freedom for the KS test are 
determined by the number of independent pairs. The 
pairs were deployed every 2 d and 100-km apart in 
the present simulation, so we can safely treat all pairs 


P ress et aL| [2007) 
tributions are so 
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Fig. 6. Relative dispersion (second moment of the pair- 
separation PDF) in summer (red) and winter (blue) for the 
NCOM pairs with various initial separations. The Lundgren 
and Richardson predictions are indicated in solid and long- 


t[ d] 


dashed curves, respectively. The short-dashed line is the long¬ 


time asymptotic Rayleigh relative dispersion. The horizontal 


Fig. 5. (top) Pair separation probability distribution func¬ 
tions (PDFs) in summer (red) and winter (blue) for various 
initial separations. The PDFs are shown at time t = 3.5 
d. Solid and dashed, and dot-dashed curves are theoretical 
PDFs in the Lundgren, Richardson, and Rayleigh regimes, re¬ 
spectively. Parameters T (with units d —1 ) and p (with units 
km 2 / 3 d — - 1 ) for the Lundgren and Richardson PDFs, respec¬ 
tively, were obtained by comparing the theoretical and ob¬ 
served dispersion curves via least squares out to the time by 
which \J ( r 2 ) = 5ro. The relative diffusivity k ,2 (with units 
km 2 d —1 ) for the Rayleigh distribution was calculated as the 
mean of the relative dispersion divided by 4 1 after t = 10 d 
(Appendix A). When a theoretical curve is statistically sim¬ 
ilar to the observed PDF in summer at the 95% confidence 
level, it is drawn in bold, (bottom) The Kolmogorov-Smirnov 
probability for the theoretical curves vs the summer PDF as 
a function of time. 


as independent. In the figure, the bold curve is sta¬ 
tistically similar to the observed summer PDF, with 
95% confidence. The lower panels show how the KS 
probability evolves with time, over the first 15 days. 
When a curve exceeds 0.05, it is plotted in bold. 

Consider the case with tq = 1 km initial separa¬ 
tion. The summer and winter PDFs are very sim¬ 
ilar and are highly kurtosed, with most pairs hav¬ 
ing small separations but some having much larger 
ones. The Lundgren distribution (solid) has a simi¬ 
lar shape, and indeed is statistically the same at the 
95% confidence level. In contrast, the Richardson 
(long-dashed) and Rayleigh (short-dashed) distribu¬ 
tions have a shorter tail and much larger mode, re¬ 
spectively. The KS probability (lower panel) though 
suggests the Lundgren PDF is significantly similar 
only during a brief period, near t = 3 d. 

The results with tq = 5 and 10 km are similar, in 
that observed PDFs resemble the Lundgren distribu¬ 
tion. The similarity moreover is significant at the 95 


dashed line indicates the deformation radius. Note that the 
t-axis is logarithmic until t = 10 d and linear thereafter. 

do not find similarity with the Rayleigh distribution 
in any case, during these initial 15 d. 

The relative dispersion curves are shown in Fig. 
[6} The winter (red) and summer (blue) curves are 
nearly identical, again supporting similar behavior in 
the two seasons. The initial growth is close to expo¬ 
nential initially, with an e-folding time on the order 
of 1 d. The exception is the tq = 30 km case, where 
the separations are only briefly below the deforma¬ 
tion radius (indicated by the horizontal line). With 
the two larger initial separations, the late dispersion 
increases linearly, in line with diffusive growth, but 
this is not apparent in the 1 and 5 km cases. There 
is a clear suggestion of Richardson-like growth with 
r 0 = 10 km, but in the other cases the dispersion is 
either greater or less than this. 

Being the fourth moment of the PDF, the kur- 
tosis (Fig. [7| is more sensitive to the tails of the 
distribution. With the smallest initial separation 
(ro = 1 km), the kurtosis grows rapidly, reaching 
values greater than 20, in both seasons. Under non¬ 
local dispersion the kurtosis grows exponentially and 
at the same rate as the dispersion (Table[l]). The ob¬ 
served growth is roughly consistent with this, with an 
e-folding time of roughly 1 d. With ro = 5 km, the 
initial growth is curtailed and the maximum values 
obtained are less. With r 0 = 10 and 30 km, the kur¬ 
tosis quickly relaxes toward 2, the asymptotic limit 
for the Rayleigh distribution. There is little support 
for a Richardson regime here; the kurtosis exceeds 
the asymptotic limit of 5.6 at the smallest separa¬ 
tions and falls below that at the larger separations. 

The velocity structure functions (S 2 ) are shown in 


% level for much of the first 10 d. The Richardson 
PDF on the other hand is not similar over the same 
period. With ro = 30 km, the PDF is statistically 
similar to the Lundgren only briefly, but then resem¬ 
bles the Richardson. However, recall that 30 km is 
only slightly less than the deformation radius, when 
the pair velocities are significantly decorrelated. We 


Fig. § for summer (solid) and winter (dashed) pairs. 
These are very similar, flattening out at scales ex¬ 
ceeding 100 km, as expected for uncorrelated mo¬ 
tion. With r 0 = 1 km, S 2 exhibits the r 2 depen¬ 
dence expected in a nonlocal regime at the smallest 
separations. With larger 7'o the curves do not have 
a unique power-law-dependence but suggest instead 
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ro = 1 km ro = 5 km ro = 10 km rn = 30 km 



![d] 


Fig. 7. Kurtosis (normalized fourth moment of the pair- 
separation PDF) in summer (red) and winter (blue) for the 
NCOM pairs with various initial separations. The Lundgren 
and Richardson curves are indicated by the solid and long 
dashed curves, respectively. Parameter-independent long-time 
asymptotic limits in the Rayleigh and Richardson regimes are 
indicated by the short- and long-dashed horizontal lines. As 
in the previous figure, the t-axis is logarithmic until t = 10 d 
and linear thereafter. 



Fig. 8. Second order structure functions in summer (solid) 
and winter (dashed). Colored curves are Lagrangian estimates 
for different initial separations. The curves with circles are Eu- 
lerian estimates. Solid and dashed straight lines have struc¬ 
ture function slopes (indicated) in the theoretical Lundgren 
and Richardson regimes, respectively. The vertical line indi¬ 
cates the Rossby deformation radius. 


a transition between the r 2 and r° asymptotic lim¬ 
its. We also calculated the Eulerian S 2 using veloc¬ 
ity time series at fixed positions (curves with circles), 
and this mirrors the Lagrangian S 2 well, supporting 
the assumptions of homogeneity and isotropy. 

Thus both the time- and distance-averaged mea¬ 
sures indicate nonlocal dispersion at the small¬ 
est separations and uncorrelated motion above 100 
km. There was some suggestion of an intermedi¬ 
ate Richardson regime (primarily in the dispersion), 
but the kurtoses and structure functions imply these 
scales are rather a transition between the small and 
large scale limits. 

One wonders of course to what extent numerical 
dissipation in the model is responsible for the nonlo¬ 
cal dispersion at the smallest scales. So we turn to 
the GLAD drifters, which are not so affected. 

4. Observed pair-separation statistics 

The GLAD pairs were obtained from quarter- 
hourly drifter positions from the GLAD experiment 




Fig. 9. As in Figs.[3](top) and[4](bottom) but based on GLAD 
trajectory pairs in three initial separation classes. 


(cf. Appendix C). The drifters were deployed near 
the DwH site, as indicated by the black dots in Fig. 
[lj Various initial separation classes were identified: 
ro ~ 0.15, 1, and 10 km. A total of 132, 127, and 
276 original pairs were obtained in each class, [j As 


Poje et al. (2014), we consider trajectory records 


spanning the initial 25 d after deployment, to avoid 
enhanced windage effects on the drifters during the 
passage of hurricane Issac. 

As the trajectories span no more than one month, 
we assume stationarity holds. Homogeneity can¬ 
not be determined with the available data, but the 
NCOM results suggest this is not an unreasonable 
assumption. Isotropy is found for separations less 
than about 100 km or so (Fig. [9j top panel). Fur¬ 
thermore, the pair velocities are correlated over the 
isotropic scales, with correlations falling below 0.5 
above 100 km (Fig. [9j bottom panel). 

The time-based statistics are shown in Fig. [To] 
The parameters for the theoretical curves were ob¬ 
tained as before, by matching the observed disper¬ 
sion from the initial separation up to the scale at 
which the rms separation \J(r 2 ) = 5ro- In principle 
one could also treat the initial separation as a free 
parameter, since a range of initial values is present, 
but we chose to set ro equal to the mean value for 
the drifters in each chosen range. 


x The specific initial separation ranges are: 0.12 km < ro < 
0.19 km, 0.84 km < ro < 1.19 km, and 9.85 km < ro < 
10.19 km. 
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Fig. 10. As in Figs. |5| (first and second rows) , [7] (third row), 
and ml (fourth row) but for separations experienced by the raw 
(red) and lowpass-filtered (blue) GLAD trajectory pairs. The 
cut-off period for the filter is 2.5 d. 



Fig. 11. GLAD trajectory pairs with initial separation ro ~ 
1 km arranged into groups exhibiting similar behavior. The 
number ( N ) of pairs in each group is indicated. 


The GLAD PDFs 3-d after deployment are shown 
in red in the top panels of Fig. 10 for the 3 initial 
separations above. The PDFs are somewhat coarser 
than for the synthetic particles as there are fewer 
pairs, but they too are peaked at the smallest sepa¬ 
rations and exhibit extended wings. 

Again, assessing the differences with the theoret¬ 
ical PDFs requires knowing the number of degrees 
of freedom, i.e., the number of independent realiza¬ 
tions. However this number was fairly small for the 
GLAD experiment, as the deployment was made over 
a limited geographical area. Drifter clusters were de¬ 
ployed very near one another, so that the distance 
between different clusters was often much less than 
the putative energy-containing eddy scale of 100 km. 
As such, many of the pair trajectories are similar. 

In fact the trajectories can be grouped into 6 
classes, as shown in Fig. 11 The numbers of pairs 


( N) in each group are indicated in the inserts, and 
these range from 13 to 29. As most of the drifters 
in each group were deployed on a single day, these 
obviously should not be considered as independent 
realizations. The exception was the N = 17 group, 
which were launched on two different days. Thus the 
number of independent realizations here is only 7, 1 
for each class and 2 for the N = 17 group. However, 
allowing for some variation in each group, we esti¬ 


mated the degrees of freedom as 3 times this, or 21. 
As such, we effectively treat each group as a triplet 
of drifters. 

With so few degrees of freedom, one cannot dis¬ 
tinguish the theoretical curves at the 95% confidence 
level. Thus both the Lundgren and Richardson PDFs 
are statistically similar over the initial 10 d, with 
all three separations. Only the Rayleigh distribution 
can be excluded during this period. 

The relative dispersion is indicated by the red 
curves in the third row of panels in Fig. [To] (the 
curves are often obscured by blue curves indicating 
the dispersion for a lowpass filtered data set, dis¬ 
cussed hereafter). The dispersion with the smallest 
initial separation is close to exponential initially but 
increases more slowly after t = 2d. With ro = 1 km 
however the dispersion follows the Richardson pre¬ 
diction over much of the first 20 d. With this sepa¬ 
ration and ro = 10 km, the dispersion asymptotes to 
linear growth at late times. 

The kurtoses are plotted in red in the bottom pan¬ 
els of Fig. |To} With ro = 0.15 km the kurtosis grows 
and then oscillates around the Richardson asymp¬ 
totic limit of 5.6. With ro = 1 km, the kurtosis 
increases more slowly, lying between the asymptotic 
limits for the Richardson and Rayleigh regimes, while 
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Fig. 12. As in Fig. pq but using raw (solid) and lowpassed 
(dashed) GLAD trajectory pairs. 


the kurtosis is close to the latter limit for most of the 
period with ro = 10 km. 

The second-order structure functions (Fig. 


12 


solid curves) on the other hand are clearly sugges¬ 
tive of a Richardson regime. With the ro = 0.15 
km pairs, the curve exhibits a power-law-dependence 
near r 2 / 3 from the smallest scales to roughly 200 km, 
consistent with the results of Poje et al. (2014). How¬ 



ever, with r o = 10 km a plateau is observed above 
100 km, in agreement with the motion being largely 
uncorrelated at those scales. 

Thus the GLAD results are ambiguous. The PDFs 
are inconclusive, due to having too few degrees of 
freedom. The dispersion curves suggest nonlocal 
dispersion at small scales while the kurtosis and 
the structure functions are more consistent with a 
Richardson regime. What causes these differences 
among the measures, and are the results actually dif¬ 
ferent than those from the synthetic pairs? 

a. Inertial oscillations 

One important difference can be seen in the tra¬ 
jectories themselves (Fig. ED-. Many of the GLAD 
drifters experience inertial oscillations: anticyclonic 
loops with a period near 1 d (the local inertial pe¬ 
riod) (cf., e.g., Gill 1982). Frequency spectra of the 


[d- 1 ] 


Fig. 13. (left) Raw (red) and lowpassed (blue) GLAD 
drifter trajectory, (right) Corresponding zonal velocity power 
spectra. 


To test this, we applied a lowpass filter to remove 
the inertial oscillations]^] The resulting spectrum and 
the corresponding trajectory are shown in blue in 
Fig. [13] The filter eliminates the peak near the iner¬ 
tial frequency while preserving the larger meandering 
motion. 

Lowpass filtering has relatively little effect on the 
time-based dispersion measures. The blue curves in 
Fig. [lO] correspond to the filtered GLAD trajecto¬ 
ries, and in all cases these mirror the results for the 
unhltered trajectories. This is because the inertial 
oscillations on nearby drifters are very similar (see 
below) and thus do not greatly alter their separa¬ 
tion. Moreover, as the drifters return approximately 
to their previous positions every inertial period, the 
net effect is small. 

The impact on the velocity structure function how¬ 
ever is greater (Fig. 


12). With the lowpass filter 


individual velocities, as in the example in Fig. |13[ ex¬ 
hibit a significant peak at the inertial period (right). 
Except for this, the energy resides primarily at the 
lowest frequencies. Consistently, the trajectory (left) 
exhibits anticyclonic loops superimposed on a larger 
scale structure. 

While inertial oscillations have a narrow frequency 
spect rum [u nless modified by background rotation 
(e.g., Kunze 1985)], the Lagrangian motion possesses 
a range of spatial scales. These are determined by the 
particle velocity, with faster-moving particles execut¬ 
ing larger loops. So the spectral profile in wavenum¬ 
ber space is broader; as such, the oscillations could 
potentially influence the separation statistics. 


(dashed curves), the energy at small scales is reduced 
and the structure function increases faster than r 2//3 . 

Inertial oscillations are only weakly captured in the 
NCOM simulation, so this effect is missed in the syn¬ 
thetic trajectories, which are furthermore computed 
from daily velocities. To test how they would have 
altered the statistics, we added inertial oscillations to 
the NCOM trajectories. This was done by modifying 
the i ^ j positions as 

Xi(t) i->- Xi(t) + A(sin uit, coswf — 1) (5) 

Xj(t) i->- X j(t) + (A + B(rij)) (sinwf, cos uit — 1). 

( 6 ) 

Here 2n/uj = 1 d, roughly equal to the local iner¬ 
tial period. The amplitude, A , was taken to be a 
random number varying over the range of observed 
loop amplitudes. The amplitude B{rij(t )) represents 
the growing difference between amplitudes on nearby 
drifters. 

The latter was chosen to mimic the behavior of 
the GLAD pairs. The inertial wave scale is generally 
much larger than the smallest pair separations here 


2 The filter was a sixth-order lowpass Butterworth filter 

with 2.5 d —1 cut-off frequency. 
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Fig. 14. (left) The root-mean-squared (rms) separation 
for highpass-filtered GLAD trajectory pairs (black) plotted 
against the rms separation of the lowpass-filtered trajectories. 
The red curve is an empirical fit used to determine the dif¬ 
ference in amplitude of the artificial near-inertial oscillations 
imposed on the NCOM trajectories. The GLAD trajectories 
belong to the ro ~ 1 km initial separation class and the NCOM 
trajectories start at GLAD positions, (right) An example of 
an NCOM pair, without (blue) and with (red) the near-inertial 
oscillations superimposed. 


(e.g., Webster|1968 

Pollard 1980 

Young and Jelloul 

1997 

Chant 

2001) 

Consistently, the high-passed 


pair velocities decorrelate on scales comparable to 
those of the full velocities (not shown), indicating 
scales of at least 100 km. As the pair velocities are 
similar for nearby drifters, the loops are correlated 
and similarly large. But as the drifters separate, the 
difference in loop amplitude grows as the velocity 
difference grows. To gauge this, we calculated the 
rms difference in the highpass filtered pair separation 
and plotted it against the rms pair separation, deter¬ 
mined from the lowpass filtered trajectories. This is 
shown in black in the left panel of Fig. [14] The rms 
difference is roughly 0.5 km initially and increases to 
a value near 2 km at separations greater than 20 km. 

Choosing B(r i:j (t)) = jj/r^(f)/r 0 with B 0 = 

3/4 yields the separation curve shown in red. This 
mimics the observed growth fairly well. The effect 
of the addition on a single pair is shown in the right 
panel, with the modified trajectories exhibiting anti- 
cyclonic loops. 

The effect on the synthetic particle statistics is 
seen in the middle panels of Fig. [15] These are for 
a set of NCOM trajectories initiated at the GLAD 
launch positions, with i' 0 = 1 km. The blue curves 
are for the unmodified trajectories and the red have 
the inertial oscillations superimposed. The PDFs, 
relative dispersion, and kurtosis are almost unaf¬ 
fected by the addition of the oscillations. But the 
structure function is significantly altered, with that 
for the modified trajectories exhibiting more energy 
at subdeformation scales. While S- 2 increases as r 2 
for the original trajectories, the dependence for the 
modified set is nearer r 2 / 3 . 



Fig. 15. Separation PDFs (top), relative dispersion (sec¬ 
ond row), kurtosis (third row), and velocity structure function 
(bottom) for: NCOM pairs at ro = 1 km separation, as shown 
in Figs.[5]p| for comparison (left); NCOM trajectories started 
from the GLAD positions with ro ~ 1 km (middle); and ac¬ 
tual GLAD trajectories with ro ~ 1 km (right). In blue are 
results based on NCOM trajectories or lowpassed GLAD tra¬ 
jectories. In red are results based on NCOM trajectories with 
near-inertial oscillations superimposed or raw GLAD trajec¬ 
tories. In green are results based on highpassed trajectories. 


Shown in green are the statistics for the inertial 
oscillations alone, obtained from the trajectories fol¬ 
lowing highpass filtering (with a cut-off of 2.5 d). 
These show that the contributions to the dispersion 
and kurtosis from the inertial oscillations are much 
less than from the low frequency components. But 
below separations of roughly 20 km, the inertial com¬ 
ponent dominates the second order structure func¬ 
tion. 

We find that adding the oscillations affects other 
distance-based measures as well, such as the rela¬ 
tive diffusivity and the FSLE (not shown). But as 
they do not significantly alter the time-based metrics, 
they cannot explain the other differences between the 
GLAD and NCOM statistics. 

b. Sampling 

The second factor influencing the GLAD results is 
the sampling strategy alluded to earlier. Most of the 
GLAD pairs were deployed very near one another 
and behaved similarly. This reduced the effective 
degrees of freedom, and the mesoscale dispersion was 
not captured well. 

The effect can be assessed by comparing in Fig.[T5| 
the statistics for the r 0 = 1 km NCOM pairs deployed 
at the GLAD locations (middle) with those from 
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the full NCOM deployment (left) and the GLAD set 
(right). The statistics for the reduced NCOM set 
are strikingly similar to those from GLAD, both the 
time- and the distance-based measures. The kurtoses 
are somewhat lower with the GLAD drifters, and the 
transition scale below which the structure function 
from the highpassed velocities crosses that from the 
lowpassed velocities is smaller (10 vs 20 km). But 
the addition or removal of inertial oscillations other¬ 
wise affects the trajectories in the same ways. And 
the results from both sets differ markedly from those 
for the full NCOM set. 

This implies the differences between the GLAD 
statistics and those of the full NCOM set may be due 
to sampling rather than to dynamics. With identical 
sampling, the NCOM model gives an accurate rep¬ 
resentation of the dispersion and the structure func¬ 
tion. This is remarkable, given that the model reso¬ 
lution is only 1 km. The reason however can be in¬ 
ferred from the frequency spectrum in Fig. [13] except 
for the peak at the inertial frequency, the spectrum 
is dominated by low frequency motions. It is these 
motions, due to mesoscale features, which dominate 
the pair dispersion, and these are captured by the 
model. This is the essence of nonlocal dispersion— 
it is controlled by structures with larger scales and 
lower frequencies. 

5. Summary and concluding remarks 

We have investigated the dispersion experienced 
by simulated and observed pairs initiated near the 
Deepwater Horizon (DwH) site in the northern Gulf 
of Mexico (GoM). The simulated separations were 
produced using synthetic pairs advected by sur¬ 
face velocities from a data-assimilative Navy Coastal 
Ocean Model (NCOM) simulation, with an effec¬ 
tive horizontal resolution of 1 km, during July 2013 
and February 2014. The observed separations come 
from drifter pairs from the Grand LAgrangian De¬ 
ployment (GLAD), conducted in July 2012. The 
investigation involved various statistical descriptors, 
namely, the probability distribution function (PDF) 
of pair separations, its second moment (relative 
dispersion), its fourth moment (kurtosis), and the 
(second-order) velocity structure function. 

The measures are consistent for the NCOM pairs, 
suggesting nonlocal dispersion at the smallest sepa¬ 
rations and diffusive dispersion at separations greater 
than 100 km, where the pair velocities are decorre- 
lated. Similar results were obtained in both winter 
and summer seasons. The results were more ambigu¬ 
ous with the GLAD pairs. While the pair motion 
is also uncorrelated at scales exceeding 100 km, the 
dispersion regime cannot be distinguished from the 


PDFs, dispersion or kurtosis. The structure func¬ 
tions on the other hand indicate Richardson disper¬ 
sion from the smallest scales to beyond the decorre¬ 
lation scale. 

Two effects impact the GLAD results. With 
such high temporal resolution (10-min sampling), the 
drifters resolve inertial oscillations. These energetic 
motions however only weakly affect pair dispersion. 
The size of the loops depend on the drifter velocity 
(e.g., Gill 1982), and since nearby drifters have nearly 
the same velocity, the size of the loops is nearly the 
same initially. As the pair separates however, the ve¬ 
locity difference grows and the amplitude difference 
increases. However, as the drifters return approxi¬ 
mately to their starting position after an inertial pe¬ 
riod, the net effect on the separation during that time 
is small. As such, the effect is greatly overwhelmed 
by mesoscale stirring. 

The oscillations do however alter distance- 
averaged measures. The velocity structure function 
is one such measure, and this is dominated by inertial 
oscillations below roughly 10 km. Filtering the tra¬ 
jectories to remove the inertial oscillations steepens 
the structure functions without affecting the time- 
based measures. Likewise, adding inertial oscilla¬ 
tions to the trajectories from the model, which has 
weak inertial variability, causes the structure func¬ 
tion to shallow, lending the appearance of local dis¬ 
persion. 

The second effect concerns the sampling in GLAD. 
As the goal was to resolve submesoscale dispersion in 
the region, the drifters were deployed in several tight 
clusters. As these spanned scales much less than the 
correlation length scale (100 km), the pairs behaved 
similarly. We found the drifters could be separated 
into 6 distinct classes, each displaying a characteris¬ 
tic path. This reduced the degrees of freedom and 
led effectively to an undersampling of the mesoscale 
stirring. Using synthetic particles deployed at the 
same locations yielded nearly identical, and equally 
ambiguous, dispersion statistics. 

The conclusion is that the 1-km model likely cap¬ 
tures the dispersion in the GLAD experiment over 
the sampled scales. This argues in favor of non¬ 
local dispersion, because the stirring is dominated 
by large scale, low frequency motions which are 
well-resolved by the model. It also supports using 
altimeter-derived geostrophic velocities to study dis¬ 
persion here, as the dominant eddies are marginally 
resolved by altimetry. Olascoaga et al. (20131 sug¬ 
gested the mesoscale circulation dominates in shap¬ 
ing the patterns formed by drifters in the GLAD ex¬ 
periment, and the present results are consistent with 
this. 
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Previously, 


LaCasce and Ohlmann (J2003) observed 

from 


exponential dispersion among “chance pairs’ 
the Surface-CUrrent and Lagrangian drifter Program 
(SCULP), from separations of 1 km up to L D (~ 50 
km). The SCULP pairs exhibited an e-folding time 
on the order of 1 d, similar to the time scales inferred 
here, and exhibited large kurtoses. LaCasce (2010) 
found moreover that the SCULP PDFs resemble the 
Lundgren distribution. Note the SCULP drifters had 
daily positions an d so were essentially de void of in¬ 
ertial oscillations. LaCasce and Ohlmann (2003) did 
not observe diffusive dispersion at super-deformation 
scales but something closer to ballistic growth, with 
the dispersion increasing as t 2 ' 2 . Given that the pair 
motion is uncorrelated above Ld, such growth most 
likely reflects shear dispersion, due to a large scale 
flow. Indeed, many of the SCULP drifters were ad- 
vected by boundary currents. 


Poje et al. (2014) presented the first analysis of 
the GLAD pair trajectories and concluded the dis¬ 
persion was consistent with the Richardson regime, 
from the smallest sampled scale (0.1 km) to several 
hundred kilometers. However, their conclusions were 
based solely on distance-based measures (the second 
order structure function and the relative diffusivity) 
which are affected by inertial oscillations. It should 
be emphasized too that a turbulence framework can¬ 
not be applied to interpret results at separations of 
hundreds of kilometers, as the pair velocities are un¬ 
correlated. 


Jullien et al. (1999) and Jullien (2003) calculated 


separation PDFs from pairs of particles deployed in 
2D turbulent flows in the laboratory. They suggested 
that the separation PDF could be fit with an empir¬ 
ical function of the form: 


P(r, t ) 



(7) 


where a = (r 2 ) 1 / 2 , and a and b are constants. They 
claimed that the same PDF applied for both the en¬ 
ergy and enstrophy cascade ranges, with slightly dif¬ 
ferent values of a and fo[^] It is straightforward to 
show that the kurtosis for this empirical PDF ([7]) is: 


(r 4 ) 

(r 2 ) 2 


9’ 

SP= 25 - 2 - 


( 8 ) 


There is no indication of such a large asymptotic 
limit in either the simulated trajectories or the 
GLAD data. So we can most likely rule out this 
type of dispersion. 

The present results serve as a cautionary note on 
using relative dispersion to deduce kinetic energy 


’’The PDF they proposed was not properly normalized. Do¬ 
ing so yields a = b 2 / 2. 


spectra. The inertial oscillations contribute to the 
spectra but do not greatly impact dispersion. So 
finding exponential relative dispersion does not nec¬ 
essarily imply steep spectra. Conversely, having shal¬ 
lower spectra at small scales, as in the atmosphere 


( Nastro m and Gage 1985) and ocean (Callies and 
Ferrari|2013), does not rule out nonlocal pair disper¬ 


sion. 

It is also worth noting that the submesoscale spec¬ 
trum here exhibits a scale separation. The La¬ 
grangian frequency spectra has peaks at low frequen¬ 
cies and at the inertial frequency, with little energy 
in between. So it is possible to partition the flow into 
mesoscale turbulence and inertial oscillations. If in¬ 
stead an inverse energy cascade were occurring from 
the submesoscales, the intermediate scales and fre¬ 
quencies would be energetic. 

The results also have implications for the design of 
dispersion experiments. Care should be taken to en¬ 
sure sufficient sampling at scales exceeding those of 
the energy-containing eddies. Otherwise one might 
obtain many similar pair trajectories, with a corre¬ 
sponding loss of statistical confidence. 
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APPENDIX A 

PDF solutions 

Solutions to the Fokker-Planck equation (|T|) have 
been derived for the turbulent inertial ranges. These 
assume that all pairs have the same initial separation, 
so that p(r, 0) = (27rr) _1 i5(r — r 0 ). Note that p is 
normalized, i.e., (r°) = 1. As noted, the solutions 
can be obtained via the Laplace transform. 
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A scale-independent diffusivity « 2 = const occurs 
when the pair motion is uncorrelated. One- and two- 
particle statistics coincide in such a case, which is 
consistent with S 2 = const (indeed, {(vi — Vj ) 2 ) = 
2(u?), which does not depend on scale). The solution 
to ( 111 ) is given by 


p(r -' ,= tob / “(S |exp,_ 


rl + r 2 
4 K 2 t 


(Al) 


p{r,t) 


1 


47TK2 1 


exp - 


4k 2 A 


(A2) 


(r 2 ) ~ 4 k 2 A, 


(A3) 


Hi 

7*2^2 


(A4) 


p{r,t) = 



x exp — 


(A5) 


p(r, t) 


n(f3t) : 


exp 


9r 2/3\ 

' 4/3 i ) 


(A6) 


which is the 2D analogue of Richardson s 1926 


lution. The relative dispersion associated with (|A5|) 
is: 


so- 


5! /4 /3tY , , /_ „ 9r, 


(O ^ Hr 


M 6,3, 



exp 


9r, 


2/3' 


4/3f 


where M(, , ) is the Rummer’s function (Graff et al. 
2015)); its long-time asymptotic limit is given by 


(A7) 


w here Iq( ) is a zeroth-order modified Bessel function 
( LaCasce|2 010). In the long-time (t n^ 1 r, r r 0 ) 
asymptotic limit, 


(r 2 ) ~ 5.2675/3 d t 3 . 


The kurtosis of (A5) is 


1 < 


which is proportional to the Rayleigh PDF. The sec¬ 
ond (raw) moment (or relative dispersion) of (A2), 


Hi 

; r 2)2 


< 5.6, 


(AS) 


(A9) 


(with the equality holding initially), while that of its 
long-time asymptotic limit (A61 is 


as expected for a normal diffusive process. The 
fourth moment normalized by the relative dispersion 
(or kurtosis) of (A2), 


r2)2 


5.6, 


(A10) 


reflecting the self-similarity of the Rayleigh PDF. 

Pair motion is correlated in the turbulent inertial 
ranges. The Richardson regime corresponds to the 
energy inertial range (both in 3D and 2D), and the 
correlated motion sustains local dispersion. With 
E oc A; -5 / 3 , the diffusivity has the form k 2 = /3r 4//3 
(Richardson 1926 Obhukov 1941 Batchelo~r||1950[ ), 
with the constant /? is proportional to the third root 
of the energy dissipation rate. The second-order 
structure function, which is the inverse Fourier trans¬ 
form of the kinetic energy spectrum, is S 2 oc r 2 / 3 
(Kolmogorov 1941). 

The solution to ([!]) is: 


which reflects the self-similarity of the Richardson 
PDF. 

Finally, in the enstrophy cascade inertial range, 
with E oc At 3 , the diffusivity is k 2 = T _1 r 2 , where 
T is proportional to the inverse cubic root of the en¬ 
strophy dissipation rate (Lin 1972). The correspond- 
ing second-ord er structure function is S 2 oc r 2 (e.g., 
B ennett||l984 ). 

The solution to ([I]) is given by: 


P(r,t) = 


4tt 3/2 (A/T) 1 / 2 ^ 


(In r/r 0 + 2A/T) 2 
4 t/T 


,_, ,_, ___ (All) 

( Lundgren|| 1981 1 |Bennett||2006 LaCasce||2010[ ). The 
relative dispersion is 


/ 2 \ 2 

(r ) = r 0 exp —, 


while the kurtosis is 


4 ) 8 1 


(r 2 ) 


.2 \ 2 = ex P 


(A12) 


(A13) 


where / 2 ( ) is a second-order modified Bessel function 
( LaCasce|2010| ). In the long-time (t /§> ft _1 r 2 / 3 , r >• 
7"o) asymptotic limit 


The Lundgren PDF (All) is lognormal and thus not 


self-similar: it gets more peaked in time, possessing 
increasingly long tails (at large scales). Note that 
the same PDF and exponential growth occurs with 
a kinetic energy spectral slope with a > 3 (Bennett 
1984 Babiano et al.||T990 ). 
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APPENDIX B 


The NCOM simulation 

Configured for the GoM, the NCOM simulation 
employs assimilation and nowcast analyses from 
NCODA (Navy Coupled Ocean Data Assimilation) 
(Cummings 20051. Forecasts are generated by sys¬ 


tems linking NCODA with regional implementa¬ 
tions (Rowley and Mask |2014 ) of NCOM (Barron 


et al. 2006). The model has 1-km horizontal reso¬ 


lution and was initiated on 15 May 2012 from the 
then operational global ocean model Global Ocean 
Forecast System (GOFS) 2.6 (Barron et al. 2007). 
Daily boundary conditions are received from the cur¬ 
rent operational GOFS using the HYbrid Coordi¬ 
nate Ocean Model (HYCOM) ( Metzger et al.||2009 ). 
The vertical grid is comprised of 49 total levels; 34 
terrain-following er-levels above 550 m and 15 lower 
z-levels. The cr-coordinate structure has higher reso¬ 
lution near the surface with the surface layer hav¬ 
ing 0.5-m thickness. The simulation uses atmo¬ 
spheric forcing at the sea surface from COAMPS 
(Coupled Ocean/Atmosphere Mesoscale Prediction 
System) (NRL 1997) to generate forecasts of ocean 
state out to 72 h in 3-h increments. The observa¬ 
tional data assimilated in these studies is provided 
by NAVOCEANO (Naval Oceanographic Office) and 
introduced into NCODA via its ocean data quality 
control process. Observations are three-dimensional 
variational (3D-Var) assimilated (Smith et al. 2011) 
in a 24-h update cycle with the first guess from the 
prior day NCOM forecast. 


APPENDIX C 


over the DeSoto Canyon area, the deployments were 
carried out in triplets, with the drifters in each triplet 
separated roughly 100 m from each other. The main 
deployments consisted of 2 clusters of 30 triplets ar¬ 
ranged in S-shaped configurations. One cluster was 
released on 20 July 2012 centered at (88.1°, 28.8°N) 
and the other cluster on 20 July 2012 at (87.6°, 
29.2°N). Each S-track spanned an area of approx¬ 
imately 8-krn x 10-km and consisted of 10 nodes 
spaced 2- to 4-km apart. Each node was made up 
of 3 equilateral triangles with 500-m side. Another 
cluster of 10 triplets arranged in a triangular configu¬ 
ration spanning an area similar to that spanned by S- 
shaped configurations was launched on 29 July 2012 
near (87.5°, 29.0°N). Two additional clusters with 20 
triplets in total were released over 30-31 July 2012 
near (89.2°, 27.8°N) inside a cyclonic eddy feature of 
about 50 km in diameter. 
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